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Double J/lp production has been observed by the NA3 collaboration in ttN and pN collisions 
with a cross section of the order of 20-30 pb. The iptp pairs measured in n~ nucleus interactions at 
150 and 280 GeV/c are observed to carry an anomalously large fraction of the projectile momentum 
in the laboratory frame, x-^^ > 0.6 at 150 GeV/c and > 0.4 at 280 GeV/c. We postulate that these 
forward t[>tjj pairs are created by the materialization of Fock states in the projectile containing two 
pairs of intrinsic cc quarks. We calculate the overlap of the charmonium states with the \udcccc) 
Fock state as described by the intrinsic charm model and find that the tt~ N — > ipi/j longitudinal 
momentum and invariant mass distributions are both well reproduced. We also discuss double J/ip 
production in pN interactions and the implications for other heavy quarkonium production channels 
in QCD. 
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Introduction 



It is quite rare for two charmonium states to be produced in the same hadronic collision. However, the NA3 
collaboration has measured a double J/tp production rate significantly above background in multi-muon events with 
7r~ beams at laboratory momentum 150 and 280 GeV/c Q and a 400 GeV/c proton beam ||]. The integrated 
ir~ N — » ijjipX production cross section, a^, is 18 ± 8 pb at 150 GeV/c and 30 ± 10 pb at 280 GeV/c, and the 
pN — > ip^X cross section is 27 ± 10 pb. The relative double to single rate, cr^,^ /o>, is (3±1) x 10~ 4 for pion-induced 
production where is the integrated single ip production cross section. 

A particularly surprising feature of the NA3 tt~ N — > ipipX events is that the laboratory fraction of the projectile 
momentum carried by the ipijj pair is always very large, x^ > 0.6 at 150 GeV/c and x^ > 0.4 at 280 GeV/c. In 
some events, nearly all of the projectile momentum is carried by the ipnp system. In contrast, perturbative gg and qq 
fusion processes are expected to produce central iptp pairs, centered around the mean value, {x^) ~ 0.4-0.5, in the 
laboratory 

The average invariant mass of the pair, (M^) = 7.4 GeV, is well above the 2m^ threshold. In fact, all the 
events have > 6.7 GeV. The average transverse momentum of the pair is quite small, pT.Tpip = 0.9 ± 0.1 GeV, 

suggesting that ipip pair production is highly correlated The proton events have a somewhat lower invariant 

mass, (M^) « 6.8 GeV. The x^ distribution for the pN — s- tptpX events has not been reported. 

There have been attempts to explain the NA3 data within conventional leading- twist QCD. Charmonium pairs can 
be produced by a variety of QCD processes including BB production and decay, BB — > ipipX H and 0(a^) ipip 
production via gg fusion and qq annihilation Li and Liu have also considered the possibility that a 2 ++ cccc 

resonance is produced which then decays into correlated ipip pairs j?| . All of these models predict centrally produced 
tpip pairs |S-n|, in contradiction to the ir~ data. In addition, the predicted magnitude of is too small by a factor 
of 3-5. If these models are updated using recent branching ratios and current scale-dependent parton distributions, 
the predicted leading twist cross sections are further reduced J8|, suggesting that an additional mechanism is needed 
to produce fast ipip pairs. 

We also note that strong discrepancies between conventional QCD predictions and experiment have recently been 
observed for rp, ip' , and T production at large px in high energy pp collisions at the Tevatron Braaten and Fleming 
pot have suggested that this surplus of charmonium production is due to the enhanced fragmentation of gluon jets 
coupling to the octet cc components in higher Fock states \ccgg) of the charmonium wavefunction. However, this 
explanation does not explain the NA3 tpip fixed target data since the color-octet mechanism still emphasizes central 
production. 

In this paper we shall explore an alternative higher-twist production mechanism in QCD where forward single and 
double ip's are created through the materialization of intrinsic cc Fock components of the pion or proton projectile 
jnj. In such states the heavy constituents tend to carry the majority of the projectile momentum since this minimizes 
the off-shell energy of the wavefunction. We find that the iprp momentum and mass distributions are reasonably well 
described by the overlap of the charmonium states with the \udcccc) Fock configuration postulated in the intrinsic 
heavy quark model. 

The QCD wavefunction of a hadron can be represented as a superposition of quark and gluon Fock states. For 
example, at fixed light-cone time, r = t + z/c, the 7r~ wavefunction can be expanded as a sum over the complete basis 
of free quark and gluon states: = J2 m m ) ^m/v- {xi 7 kT,i,Xi) where the color-singlet states, |m), represent the 

Fock components \ud), \udg), \udQQ), etc. The boost-invariant light-cone wavefunctions, i/) m / n - (xi, kr,i, A;), needed 
to compute probability distributions, are functions of the relative momentum coordinates Xi = kf~ / ' P + and &t,i- 
Momentum conservation demands X)"=i Xi = 1 and ^T,» = 0, where n is the number of partons in a Fock state 

|m). When the projectile scatters in the target, the coherence of the Fock components is broken, its fluctuations can 
hadronize, forming new hadronic systems from the fluctuations []l2| . For example, intrinsic cc fluctuations can be 
liberated provided the system is probed during the characteristic time, At = 2pi a b/M?-, that such fluctuations exist. 

Microscopically, the intrinsic heavy quark Fock component in the 7r~ wavefunction, \udQQ), is generated by virtual 
interactions such as gg — * QQ where the gluons couple to two or more projectile valence quarks. The probability 
for QQ fluctuations to exist in a light hadron thus scales as ^(toq) / 'uiq relative to leading- twist production JL3| . 
Therefore, this contribution is higher twist, suppressed by 0(1/itiq) compared to sea quark contributions generated 
by gluon splitting. In (1 + 1) QCD the Fock state representation of each hadron can be computed explicitly, including 
its intrinsic QQ configurations, by diagonalizing a discrete form of the light-cone Hamiltonian 14 1. 

In general, the dominant Fock state configurations are not far off shell and thus have minimal invariant mass, 

M 2 = X)j to t il ' x i where tot,« = i + mf is the transverse mass of the i th particle in the configuration. Intrinsic 

QQ Fock components with minimum invariant mass correspond to configurations with equal rapidity constituents. 
Thus, unlike sea quarks generated from a single parton, intrinsic heavy quarks tend to carry a larger fraction of 
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the parent momentum than the light quarks JTl| ]. In fact, if the intrinsic QQ coalesces into a quarkonium state, 
the momentum of the two heavy quarks is combined so that the quarkonium state will carry a significant fraction 
of the projectile momentum. It was shown that large xp virtual cc or lepton pairs can be liberated by a relatively 
soft interaction with a light quark component of the projectile fl^| . For soft interactions at momentum scale /i, the 
intrinsic heavy quark cross section is suppressed by a resolving factor oc /i 2 /tv/q ]l3| - 

There is substantial circumstantial evidence for the existence of intrinsic cc states in light hadrons. For example, 
the charm structure function of the proton measured by EMC is significantly larger than predicted by photon-gluon 
fusion at large xsj Jl5f|. Leading charm production in ttN and hyperon-TV collisions also requires a charm source 
beyond leading twist |l3l , |il| . The NA3 experiment has also shown that the single J /if; cross section at large xf is 
greater than expected from gg and qq productio n ]l7t . Additionally, intrinsic charm may account for the anomalous 
longitudinal polarization of the J/ip at large Xp |18[| seen in nN — ► J/ipX interactions. 

Over a sufficiently short time, the pion can contain Fock states of arbitrary complexity. For example, two intrinsic 
cc pairs may appear simultaneously in the quantum fluctuations of the projectile wavefunction and then, freed in 
an energetic interaction, coalesce to form a pair of i/j's. We shall estimate the creation probability of \nvcccc) Fock 
states, where ny = du for tt~ and ny = uud for proton projectiles, assuming that all of the double J/tp events arise 
from these configurations. We then examine the x^ and invariant mass distributions of the ipip pairs and the x<p 
distribution for the single ip's arising from these Fock states. 



Intrinsic Charm Fock States 

The probability distribution for a general n-particle intrinsic cc Fock state as a function of x and kr is written as 



nr=i dxid?k Tit " s cc (to 2 - Yli=i( m T,i/ x i)) 2 

where N n normalizes the Fock state probability. In the model, the vertex function in the intrinsic charm wavefunction 
is assumed to be relatively slowly varying; the particle distributions are then controlled by the light-cone energy 
denominator and phase space. This form for the higher Fock wavefunctions generalizes for an arbitrary number of 
light and heavy quark components. The Fock states containing charmed quarks can be materialized by a soft collision 
in the target which brings the state on shell. The distribution of produced open and hidden charm states will reflect 
the underlying shape of the Fock state wavefunction. 

The invariant mass of a cc pair, M C c, from such a Fock state is 

dP f " dx - dP 

m% = J n **** nr =1 ^W, 5{Xc ~ c ~ Xc ~ Xz) (2) 

x 5(kr,c + k T< c - k T , cZ ) S 
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where n = 4 and 5 is the number of partons in the lowest lying meson and baryon intrinsic cc Fock states. The 
probability to produce a J/ijj from an intrinsic cc state is proportional to the fraction of intrinsic cc production below 
the DD threshold. The fraction of cc pairs with 2to c < M C c < 2to,£> is 



c/h 



im D r ip. / r s rip. 

dM cc S" / / dM cc 4tt ■ ( 3 ) 
4m 2 dM% / J im 2 dM% 



The ratio f C c/ir is approximately 15% larger than f C c/ p for 1.2 < m c < 1.8 GeV. However, not all cc's produced below 
the DD threshold will produce a final-state J/tp- We include two suppression factors to estimate J/ip production, 
one reflecting the number of quarkonium channels available with M C c < 2mo and one for the c and c to coalesce with 
each other rather than combine with valence quarks to produce open charm states. The "channel" suppression factor, 
s c sa 0.3, is estimated from direct and indirect J/tp production, including xi and Xi radiative and ip' hadronic decays. 
The combinatoric "flavor" suppression factor, s/, is 1/2 for a \udcc) state and 1/4 for a \uudcc) state. In Fig. 1 we 
show the predicted fraction of ^'s produced from intrinsic cc pairs, 

fi/j/h = s cSff C c/h , (4) 
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as a function of m c . We take m c = 1.5 GeV, suggesting ss 0.03 and f^/ p ~ 0.014. 

The nuclear dependence arising from the manifestation of intrinsic charm is expected to be a a ~ &nA 2 / 3 , charac- 
teristic of soft interactions. Indeed, in their single J/ip measurements [jF7j, the NA3 collaboration defined a separate 
"diffractive" cross section at large xp scaling like A - 77 for it A and A"^ for pA interactions. We thus identify the 
diffractive cross section with af c , as in our previous study p9| . The intrinsic cc production cross section from an 
\nvcc) configuration, 

a ic (hN)=P ic a%,^, (5) 



was estimated in Ref. Jl3|] to be <Ji c (irN) « 0.5 fib and a- lc (pN) « 0.7 fib for a beam momentum of 200 GeV. At low 
xf, a different, nearly linear, A dependence was observed. The ratio of diffractive to total J/ifi production can then 
be determined from the NA3 data jl?]]. The soft interaction scale parameter fi 2 ~ 0.2 GeV 2 , is thus fixed by the 
assumption that the diffractive fraction of the total production cross section is the same for charmonium and charmed 
hadrons. 

The diffractive J/ip cross section is (18 ± 3)% of the total it A — > J/ip production cross section and (29 ± 6)% of 
the pA — > J/ip production Jl?]]. The cross sections in the J/ijj — » fi + fi~ channel are Bo^(tx~ N) = 6.5 ± 0.6 nb at 150 
GeV and Ba^{pN) = 3.6 ±0.6 nb at 200 GeV. Removing the branching ratio 0, then implies af c (-K~N) = 15.7 ±3.0 
nb and af c (pN) = 14.3 ± 3.3 nb. The proton cross section should be regarded as an upper bound since the bulk of 
this cross section is at xf < 0.2 where the intrinsic charm contribution is expected to be small. Using the data for 
xp > 0.2 suggests a lower bound, af c (pN) = 5.6 ± 1.3 nb. We relate af c to the intrinsic cc production cross section by 

4(hN) = U /h <r ic (hN) . (6) 

Using the estimates of (Ti C (7rV) from Ref. |l3| with the values of f^/h in Eq. (4), we find af c (ir~ N) — 15 nb and 
af c (pN) = 9.8 nb, in agreement with the NA3 single ip data]]. 



Double J/tp Production from Intrinsic Charm Fock States 

If one assumes that all of the NA3 double J ftp events arise from intrinsic \nvcccc) Fock states, then the required 
normalization for this state can be determined from 

<rf(hN) = fl /h ^ <? ic (hN) = U /h at(hN) (7) 

Mc MC 

where P cc is the probability to produce a pair of intrinsic cc states in the projectile. Then = afj' (n~ N) « 20 pb 

0, requiring p cc w 4.4% p c . If p cc is independent of the projectile, then we predict afJ'ipN) w 6 pb, which is 20% 
of the measured cross section H, but consistent with our previous calculations (ll|. If, instead, the ratio a^/a^ Q| 
is assumed to be independent of the projectile, the cross section increases to 14.6 pb, more compatible with the data 
and Pi cc « 10.6% P; c . Finally, if we assume afJ'(pN) = 27 pb, as suggested by the data of Ref. Jl[, then we require 
p cc ~ 20% P; c . Note that the probability of making a second intrinsic cc pair is relatively easy once one cc is present 
in the Fock state. 

The Tptp invariant mass distributions measured by NA3 [0^] are shown in the histograms of Fig. 2. The pair mass 
distribution predicted from the \nvcccc) Fock state is 

J§r = / n n ^ d<^ 3 m j^ kTi (8) 

x S — — ^ I S(k T . Cj + fc Tj c 3 - kT^ ] )8(x^ j - x Cj - x-c 3 ) 




r / ^TM m T,i: 1 '•'•lew i 
X d ] 0{k T ,4, 1 + Kt,i/j 2 - kT,inp)0(X^ - X t f >1 - x^ 2 ) 

\ x^pip x^ 1 



1 We use B(J/ip — > fi + fi ) ~ 7.4% for consistency with the NA3 measurements. If B(J/ip — > /i + fi ) ~ 6% is used instead, the 
single cross sections increase by 20%. 
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where n = ny + 4 and 



dPiCC = N n [ai(M c ,)f ^V^-^h^ (9) 



n™=i dxid 2 k T ,i " s Cc (m 2 h - T,7=i( m T,i/ x i)¥ 

The delta functions insure conservation of momentum for both J/ip mesons and the ipip pair. Here, 2m c < < 2tod 
and m c = 1.5 GeV. In Fig. 2 the predictions for the ipip pair mass distributions are compared with the tt~ N (a) and 
pN (b) data and normalized to the data. Without any kx dependence, the mass distribution is strongly peaked 
at threshold. The kx dependence smears out the pair distribution, increasing (M^) by 0.5 — 1.0 GeV. We find 
(M^) « 7.7 GeV for the pion beam and (M^) w 7.4 GeV for the proton beam. The intrinsic charm model is very 
successful in reproducing the strongly correlated features of the data. The smaller value of (M^) in pN interactions 
also agrees with the trend of the data [||| . 

Although the kx dependence is needed to calculate the mass distributions, it is sufficient to use a mean value 
of k T to calculate the x distributions since the shapes of the longitudinal momentum distributions are essentially 
independent of the transverse momentum. Then the x distributions from an \nvcccc) state, following Eq. (1), can be 
written as 

dPicc - N n [a*(Mcc)? , /^H^^ > (10) 



dxi---dx n " s cc (ml - Y^=i( m "i I Vi)) 2 



where again n = ny + 4 and m t = y mf + (fc^J is the average transverse mass. Assuming {k T t ) is proportional to 
the square of the constituent quark mass, we adopt the effective values fh c = 1.8 GeV and fh q = 0.45 GeV, as in our 
previous work |l9|,Q. The distribution for a single J/ip in a \nvcccc) state is 

= [ T\ dxidx^ 2 dPicc — ^ _ x _ x _ )<J(aty - x C2 - Xc 2 ) ■ (11) 
«^Vi J i dxi...dx n ' 

We find (x^p) = 0.36 for the pion projectile and (a;^,) = 0.33 for the proton. The number of single J/tft's is twice the 
number of ipil) pairs. The single J/ip distributions from the \nvcccc) state have a lower average x^ than those from 
\nvcc) Fock states, where (x^) = 0.62 for a pion and (x,^) = 0.51 for a proton The pair distribution can be 

computed from 

A D r dPi cc 

dx^ 1 5(x^ - x^ 1 - Xfa) . (12) 



dx^ptp J dxip 1 

The intrinsic charm model predicts (xm) — 0.72 for the pion and (x^) — 0.64 for the proton. We compare the 
frame-independent calculation to histograms of the combined 150 and 280 GeV data assuming x^ = P^ b /pbeam in 
Figure 3. The tptp pair distributions are shown in Fig. 3(a) and 3(c) and the associated the single J/ip distributions 
in pair events are shown in Fig. 3(b) and 3(d). Both are normalized to the data with the single J/ip normalization 
twice that of the pair. 



Other Tests of the Intrinsic Heavy Quark Mechanism 

The intrinsic charm model provides a natural explanation of double J/ip hadroproduction and thus gives strong 
phcnomcnological support for the presence of intrinsic heavy quark states in hadrons. While the general agreement 
with the intrinsic charm model is quite good, the excess events at medium x^ suggests that intrinsic charm may not be 
the only ipip QCD production mechanism present or that the model parameterization with a constant vertex function 
is too oversimplified. The x^ distributions can also be affected by the A dependence. Additional mechanisms, 
including an update of previous models , will be presented in a separate paper Q . 

The intrinsic heavy quark model can also be used to predict the features of heavier quarkonium hadroproduction, 
such as TT, Tip, and (bc)(cb) pairs. Using m& = 4.6 GeV, we find that the single T and TT pair x distributions 
are similar to the equivalent ipip distributions. The average mass, (Mxt), is 21.4 GeV for pion projectiles and 21.7 
GeV for a proton, a few GeV above threshold, 2mx = 18.9 GeV. The x~f%p pair distributions are also similar to the 
ipip distributions but we note that (xt) = 0.44 and (x^) — 0.30 from a \udccbb) configuration and (xt) = 0.39 and 
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(xj/,) = 0.27 from a \uudccbb) configuration. Here (My^) = 14.9 GeV with a pion projectile and 15.2 GeV with a 
proton, again a few GeV above threshold, mx + = 12.6 GeV. 

It is clearly important for the double J/tp measurements to be repeated with higher statistics and also at higher 
energies. The same intrinsic Fock states will also lead to the production of multi-charmed baryons in the proton 
fragmentation region. It is also interesting to study the correlations of the heavy quarkonium pairs to search for 
possible new four-quark bound states and final state interactions generated by multiple gluon exchange 0. It has 
been suggested that such QCD Van der Waals interactions could be anomalously strong at low relative rapidity p2| , p3| . 

There are many ways in which the intrinsic heavy quark content of light hadrons can be tested. More measurements 
of the charm and bottom structure functions at large x f are needed to confirm the EMC data Q . Charm production 
in the proton fragmentation region in deep inelastic lepton-proton scattering is sensitive to the hidden charm in the 
proton wavefunction. The presence of intrinsic heavy quarks in the hadron wavcfunction also enhances heavy flavor 
production in hadronic interactions near threshold. More generally, the intrinsic heavy quark model leads to enhanced 
open and hidden heavy quark production and leading particle correlations at high in hadron collisions with a 
distinctive strongly-shadowed nuclear dependence characteristic of soft hadronic collisions. 
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Figure Captions 



Figure 1. The fraction of intrinsic cc states below the DD threshold that produce J ftp's is given as a function of 
charmed quark mass for pion (solid) and proton (dashed) projectiles. 

Figure 2. The calculated tptp mass distributions, Eq. (8), are compared with histograms of the NA3 data jl],[|] from 
7r~ N (a) and pN (b) production. The calculations are normalized to the total number of tptp events. 

Figure 3. The tptp pair distributions are shown in (a) and (c) for the pion and proton projectiles. Similarly, the 
distributions of J ftp's from the pairs are shown in (b) and (d). Our calculations are compared with the ir~N data at 
150 and 280 GeV/c Q|. The distributions are normalized to the number of pairs from both pion beams (a) and 
the number of pairs from the 400 GeV proton measurement (c). The number of single J ftp's is twice the number of 
pairs. 
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